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SUMMARY 


This volume of the report on The Gordian Refinery Simulation Model 
contains the detailed systems and prograiranlng documentation for the 
computer program. The program predicts the flow streams* material* 
energy and economic balances for a refinery processing shale oil, coal 
oil, and petroleum crudes with emphasis on the production of jet fuel of 
varying end point and hydrogen content specifications. This volume 
Includes a detailed description of all program modules, subroutine 
common blocks and variables, program function statements and data 
values, files used during program execution and a list of program limitations 
as determined by the size of dimensioned arrays. The first two volumes 
of this report and a complete Fortran listing of the program are necessary 
supplements to this volume. 
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1.0 IMTRQDUaiON 


Major prica increases and the inipending shortage of petroleum 
reserves with respect to increasing product demand has brought about a 
serious examination of possible changes in jet fuel composition. 
Specification aviation turbine fuel (AS'fM D-1655) is produced from mid- 
distillate petroleum fractions, which compete with ever growing demands 

i 

for diesel, fuel oil, and petrochemical feedstocks. Increased distillate 
production from present crudes is feasible, but conversion of gas oils 
and residuals increase the aromatic content of the mid-distillate pool. 
Moreover, promising alternate crude sources, such as shale oil, tar 
sands, and coal liquids yield distillates also with increased aromatic, 
nitrogen and sulfur contents. Special processing would be required to 
produce present specification aviation turbine fuel from these sources. 

This view of the future has stimulated a reexamination of the 
optimum combination of jet fuel specifications, with respect to the 
refinery processing, the supply distribution system, the aircraft fuel 
system and the fuel combustion qualities. The goals of current studies 
are assessing the suitability of jet fuels produced from cracked petroleum 
and alternate crude sources and developing a data base which will allow 
optimization of future fuel characteristics. Future aviation turbine 
fuel specifications must represent a trade-off between energy and cost 
efficiency of manufacture and aircraft and engine design and performance. 

This report deals with the refinery portion of the overall program. 

In order to have a systematic way of determining the energy efficiency 
of the production of various product slates involving different crude 


Z 


sources and different processing schepies, the Lewis Research Center of 
NASA has supported the development of this computer model for petroleum 
refinery operation. The primary objectives of this model are: 

1. The flexibility to configure a refinery involving any or all 
of the process units commonly employed in the production of 
gasoline, jet fuels, and mid-distillates; 

2. The ability to produce jet fuel blends of varying end-point 
specification and varying specified hydrogen content as part 
of the total slate of products; 

3. The ability to handle synthetic crudes (shale and coal derived) 
with varying severities of hydroprocessing; 

4. The determination of overall refinery energy efficiency; 

5. The determination of sulfur, nitrogen, and hydrogen material 
balances for each process unit and for the overall refinery; 
and 

6. The capabilily of carrying out economic calculations. 

The Gordian Refinery Simulation Model, presented herein, has all 
the above capabilities. This report is the third of three volumes. 

Volume I (NASA CR-1 35333) is a detailed description of the program, 
input data, and sample output; and Volume II (NASA CR-135334) contains a 
description of program data bases and correlations. The complete 
documentation and program tape are available through the Computing 
Software and Management Information Office (COSMIC) under the number 
LEW-13047. 

This volume is designed to aid the user in understanding the 
construction of the program ix> the point where a programmer could be 
able ta make program changes in order to extend the scope of the calculation, 
alter data base values, update program correlations and increase the 
accuracy of the program results. 


Exhibit 1 on the following page shows the program logic and the 
interrelationship between the main program and its subroutines. Reference 
is made in this volume to the preceding two volumes of this report and 
to the Fortran program listing. These are therefore required for a 
complete understanding of the contents. 
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EXHIBIT 1 - Block Flow Diagram of Program Loatc 
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2.0 MAIN PROGRAM AND SUBROUTINE DESCRIPTIONS 


The following sections contain the detailed descriptions of the 
main program and the twenty subroutines which comprise the Computer 
Model for Refinery Operations. Section 2.1 describes the functions of 
the main program and each of the subroutines and the nature rf the 
calculations performed within each of the program modules. Section 2.2 
describes the subroutine calling sequence and the subroutine heirarchy. 


2.1 Detailed Description of Program Modules 

An overall block flow diagram describing the logic of the Gordiar 
Model for Refinery Operations is shown in Exhibit 1. A description of 
the sequence of calculations and operations performed within each program 
module follows: 

A. MAIN PROGRAM 


The main program performs the following operations and computations 


in sequence: 

(1) Reads changes to the crude oil assay data base and the process 
unit yield and properties data base. These changes remain in 
effect for the entire run sequence and are stored on temporary 
disk files with logical designations of 16, 17 and 18, A 
detailed description of these files and their contents is 
given in Section 5.0. 

(2) Commencing with main program statement 777, the Initializing 
subroutines INITB and INIT\f are called to reinitialize variable 
values for each case. 

(3) Calls the process unit yield and properties data base sub- 
routines (DBASE) and the crude oil assay data base (CBASE). 

This brings all of the model data base values into core, 

(4) Reads the problem specific information for each case of the 
entire run sequence. Base case information is written and 
stored on disk fl'^e 20 while subject case changes to the base 
case are written on disk file 19. 
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(5) Connnencing with program statement 6, changes to the crude oil 
assay and process unit yield and property data base are read 
from disk files 16, 17 and 18 to overlay the CBASE and UBASE 
subroutines data base values with those supplied by the program 
user. These changes remain In force throughout the entire job 
step. 

(6) If a base case Is being processed the input data is next read 
from disk file 20. If a subject case 1s being processed, the 
file 20 read Is followed by a file 19 read, thus overlaying 
the base case with the subject case changes. 

(7) Subroutine CONAPI is called to convert all of the Input API 
gravity values to specific gravities. 

(8) Commencing with statement 888 the properties of shale oils, 
coal oils, and petroleum based oils Input to the refinery are 
blended Individually according to volume or weight as appropriate. 
The main program function statements BFPF and FPvFP are called 
within this calculational sequence to blend freezing point by 
blending index values. This computational sequence extends to 
the statement prior to main program statement 9, where the 
program loops back to statement 888 to do the blending calculation 
for coal oil and petroleum based refinery feeds (shale oil 
properties are blended in the first pass thru the loop). 

Subroutine INITV is called at the beginning of this loop in 
order to reinitialize various intermediate variables which are 
used in the blending calculations. 

(9) Commencing with main program statement 9, the composite volumes, 
specific gravity, sulfur, nitrogen and hydrogen contents of 

the shale oil, coal oil and petroleum based refinery feeds are 
blended in order that the overall refinery mass, sulfur, 
nitrogen and hydrogen balances may be calculated at a later 
stage. The total heat content of the composite refinery feed 
streams Is also computed for later use in the calculation of 
the overall refinery energy efficiency. 

(10) Commencing with the DO- LOOP statement 16 individual process- 
unit vol^_e;M£ balances are calculated. These process unit 
computations extended for about three hundred and seventy 
statements until statement 133. 

The locations of the calculations for each process unit are specified 
below: 


orb&inax. Page zb 
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statement 

Number 


Process 

Unit 


Seven prior to 16 

Kerosene hydrotreater (petroleum) 

200 

Shale oil kerosene hydrotreater 

208 

Coal oil kerosene hydrotreater 

13 

Coker 

22 

Fluid catalytic cracker 

33 

Thermal gas oil cracker 

34 

Gas oil hydrocracker (petroleum) 

Z.6 

Shale oil gas oil hydrocracker i 

224 

Coal oil gas oil hydrocracker 

44 

Gas oil desulfurizer (petroleum) 

55 

Distillate desulfurizer (petroleum) 
Distillate hydrocracker (petroleum) 

66 

77 

Catalytic naphtha reformer 

88 

Butane isomerizer 

99 

Alkylation unit 

111 

Polymerization 

122 

Hydrogen unit 


(n) Commencing with statement 133, the effect of hydrotreated 
kerosene sales, desulfurized diesel sales and raw gas oil 
sales on the refinery balance and the availability of fuel 
blending stocks is computed. This calculatlonaT sequence 
extends thru program statement 150. 

(12) Subroutine UBASE 1s recalled to bring into computer storage 
the specific gravity, freezing point and PNA "deltas". These 
deltas correspond to the change In properties between the 
product and feed streams for the specific process units. 

Files 16 and 18 are reread to overlay these values with those 
{if any) which are user specified for the run sequence. 

(13) Subroutine MATCLC is called to calculate the individual process 
unit mass balances, along with the sulfur, nitrogen and hydrogen 
balances. A lengthy sequence of program calculations is thus 
Initiated which take place outside of the main program. 

During this sequence of calculation, optional reports are 
generated which give the mass, sulfur, nitrogen and hydrogen 
balances for each process unit along with the energy efficiency 
report for each unit within the user specified refinery configuration. 

(14) Subroutine BLENDF Is next called. Again a lengthy chain of 
calculations 1s Initiated outside of the main program which 
calculates the volumes, masses and all properties of the jet 
fuel, distillate and residual blends produced by the oil 
refinery. Blend reports are generated during this sequence of 
calculations. 
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(15) The final segment of the main program thus has in hand all of 
the process unit material balance and fuel blend information 
necessary to calculate the overall refinery balances and to 
generate the corresponding reports. The sequence of calculation 
and operations within the final segment of main program is as 
follows: 

(al Calculate all refinery inputs 

fb) Calculate light end products 

(c) Calculate excess naphtha product 

(dj Calculate the jet fuel, product blend 

(e) Calculate the distillate fuel product blend 

ffj Calculate the residual fuel product blend 

fgj Calculate the coke sales 

(h) Calculate the kerosene sales (side stream) 

(i) Calculate the gas oil sales (side stream) 

(j) Calculate the diesel oil sales (side streams) 

(k) SuRsnarize the refinery process unit operations 

as elements of capacity and feed volume arrays 

(l) Do the refinery gasoline bland calculations 

(m) Generate the overall material balance report 
(crudes followed by products) 

(n) Generate the refinery process unit operations reports 

(o) Generate the refinery energy efficiency summary 
report 

(p) Generate the refinery fuel, power and steam summary 
report 

(q) Generate the gasoline pool report 

(16) Subroutine ECON is called at the end of the main program to do 
the optional refinery economic calculations. Following this 
set of calculations and the generation of the economic summary 
calculations, a return is made to the main program. The "GO 

TO 777" statement in the main program loops back to the beginning 
of the main program in order to begin the calculation for the 
next case within the job sequence. Note that the read statements . 
which pertain to overlaying the crude oil data base and the 
process unit yield and properties data base are not repeated 
since these are performed at the beginning of the entire run 
sequence and remain in force throughout the entire run. 

Separate jobs must be submitted if it is desired to make 
alternative data base changes. 

B, BLOCK DATA 


The BLOCK DATA subroutine initializes several arrays and variables 

prior to execution of the program. This initialization of 
course occurs at the beginning of the execution of the run sequence. 
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Various subroutines (later described) are used to reinitialize key 
variables at the beginning of the execution of successive case runs 
within a run sequence. The following arrays are initialized in BLOCK" 

DATA: 

(1) All refinery process unit yield and property arrays 

! Z) All refinery crude oil volumes and all crude oil array values 
3) The case title is initialized as blanks 
4 j All intermediate refinery stream volumes 
5) All refinery fuel bland volumes and properties 

In addition all of the refinery stream names are set in array 

SNAMES, all refinery process unit names in UNAMES and all crude oil 

names in array CN. 

Certain economic subroutine values are also initialized in BLOCK 
DATA. These include the capital investment costs of the process units 
(array UNTCPC), the BPD capacities to which these dollar investments 
correspond (array UNTSTR)» the power law coefficients -for 
investment (array UNTEXP), the process unit dollar per barrel Investment 
costs (array UNTVC), the dollar per day labor costs for each process 
unit (array UNTLC) and the annual process unit maintenance costs expressed 
as a percentage of the capital investment (array UKTWP). BLOCK DATA • 
also initializes the refinery stream day factor to 0.92 and various mass 
and energy balance accumulators are initialized to zero. 

C. INITB 

This subroutine initializes select input variables v 'a new base 
case is specified within a run sequence. The classes of variables which 
are initialized include the process unit capacities, process unit severities, 
process unit feed volumes (which are calculated during the course of 
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previous runs) jet and other fuel blend specifications (these are initialized 
and reinitialized at zero to suppress the blend production' unless specified 
in the input), fuel blend and miscellaneous by product stream productions 
(initialized at zero), various mass and energy balance accumulators 
(initialized at zero) and finally the report writer option (IREP) which 
is initialized at 1 for limited report generation, 

D. INITV(IJ) 

The purpose of this subroutine is to reinitialize accumulators 
prior to each run .along with refinery stream and blend property arrays. 

The entire subroutine is invoked by calling INITV(l). A limited number 
of variables are reinitialized by calling INITV(2) - this is done at 
various points within a case calculation to reset specific accumulators, 

£. CBASE 

This subroutine contains the crude oil assay information for the 
first nineteen crude oils. A complete explanation of the crude array 
variables is given on pages 27 thru 29 of Volume I, 

F. CBASEl 

This subroutine contains the crude oil assay information for the 
twentieth thru the twenty-sixth crude oils. The split into two subroutines 
is made in order to accomodate the compilation requirements of certain 
IBM computer operating systems. 

G. UBASEfli n 

This subroutine contains the following refinery process unit yield 
and stream quality arrays; 

(1) Process unit yields as a function of the feed stream 
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(2) Sulfur and nitrogen contents of refinery process streams - 
these are expressed as a function of the feed stream contents. 

(3) Octane values at 0 and 3 ml of tetraethyl lead for the gasoline 
range refinery streams. 

(4) Viscosity blending values of distillate and residual fuel. 

(5) The specific gravity and heat of combustion of distinct 
light ends species such as hydrogen. 

The above arrays are established within the program at their data 
base values by calling UBASE(l). The call to UBASE(2) establishes 
’’deltas” for various properties, where the deltas are the differences 
between the process unit product and feed stream properties. These 
Include changes in specific gravity, freezing point, and paraffin, aromatic 
and napthene contents. In the case of gas oil feed streams being charged 
to process units which produce distillate streams, absolute rather than 
differential freezing point and PNA properties are specified as a function 
of the process unit severity level. 

A complete explanation of all UBASE process unit yield and physical 
property variables Is given on pages 29 thru 45 of Volume I and Exhibits 
5 and 6 of that volume. 

H. MATCLC 

This subroutine converts the volumetric stream balances described 
in the Main Program to mass balances and computes the sulfur, nitrogen 
and hydrogen balances. Two function statements are referenced repeti- 
tively within MATCLC to convert from specific to API gravity (TAPI) and 
from API to specific gravity (TSP6G). Both are required since the feed 
stream specific gravities for each process unit must be converted to API 
gravity in order to add the data base "delta" API gravity between product 

and feed streams and this sum must then be reconverted to specific 

gravity in order to calculate the mass balances. 
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The Initial segment of MATCLC thru statement 141 sets the process 
unit gravities for the thirteen light end species to the correct values. 
The logic of the remainder of the subroutine is repetitive in nature, 
with successive sections calculating the mass, sulfur, nitrogen and 
hydrogen balances for each process unit specified for the oil refinery 
configuration. 

Four subroutines are called for each pipocess unit calculation: 

ISTN - sets the process unit stream numbers according to the 
stream codes listed in Exhibit 5 of Volume I. 

HYDCAL - calculates the hydrogen content of each stream as a 
function of API gravity, mean boiling point and sulfur and nitrogen 
content. The sole exceptions are for the petroleum, shale and coal oil 
crude units since the hydrogen contents for these process unit streams 
are obtained directly from the crude oil assays (subroutines CBASE and 
CBASEl) and It is hence not necessary to call HYDCAL. Subroutine HYDCAL 
returns to MATCLC with a calculated hydrogen content for each stream vs. 
the value of 0.0 that was originally transmitted from MATCLC. Program 
initialization sets these values at 0.0 unless overridden by user values. 
•These overides may be provided for selected streams involved in the 
blending of aviation turbine jet fuel as listed on pages 39 and 40 of 
Volume I. 


UHTREP - generates the process unit material balance report. 

UHTENR - generates the process unit energy consumption report. 

The specific sections of MATCLC where each process unit calculation 
is performed are listed below in the order of calculation: 


Statement Number Process Unit 


141 

80 

151 

1 

154 

81 

4 

5 

6 

157 

9 

10 


Petroleum crude unit 

Coal oil crude unit 

Shale oil crude unit 

Petroleum kerosene hydrotreater 

Fluid coker 

Visbreaker 

Fluid catalytic cracker 

Thermal cracker 

Petroleum gas oil hydrocracker 

Petroleum gas oil desulfurizer 

Petroleum middle distillate desulfurizer 

Petroleum middle distillate hydrocracker 


T5 



statement Number 


Process Unit 


6161 

Shale oil kerosene hydrotreater 

176 

Coal oil kerosene hydrotreater 

177 

Shale oil gas oil hydrocracker 

178 

Coal oil gas oil hydrocracker 

nil 

Naphtha reformer 

12 

Butane isomerizer 

13 

Alkylation unit 

14 

Polymerization unit 

15 

Hydrogen plant 


The material balances calculated ‘by subroutine MATCLC are an 
integral part of the overall program calculations. The UNTREP and 
UNTENR report subroutines which are called from MATCLC return immediately 
without issuing reports if the report option indicator IREP is not set at 
2 or higher. 

I, ISTN (NS, II, 12... 115) 

This subroutine performs an important function in terms of the 
overall program logic. The calls to this subroutine are made from the 
main program and subroutine MATCLC in order to identify the coded stream 

CD 

numbers for subsequent strings of volumetric and mass balance calculations 
The first argument NS is the total number of stream numbers to be coded, 
and n thru 115 are the coded stream numbers themselves. Thus up to 
fifteen stream numbers may be sat, corresponding to (for example) two 
process unit feed streams and up to thirteen product streams. This is 
more than sufficient for any of the process units represented in the 
model. Excess arguments are coded with a dummy variable. For example, 
the call statement for the alkylation unit mass balance from MATCLC is 
as follows: 


(1) A full list of the coded stream numbers is given in Exhibit 5 of 
Volume I. 
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Call ISTN (4,8,10,2*1 

where from Exhibit 5, Volume I, 

Stream 8 is isobutane (a reed) 

Stream 10 is butylene ( a feed) 

Stream 21 is light alkylate (a product) 

Stream 22 is a heavy alkylate (a product) 

The total number of streams involved is four, as indicated by the 
first call argument. M is a dummy variable which is defined in MATCLC. 
Feed streams are listed first to accomodate the reporting functions 
handled by subroutine UNTREP. 

Once the series of coded stream numbers is set within subroutine 

ISTN, the series of stream numbers is transmitted to the other program 

subroutines through labeled common/ ISTNUM/ which contains the array 

ISTR(66) and NSTOT. In the example of the alkylation process unit, the 

XSTR array positions will contain the following: 

ISTR(l) 8 
ISTR(2) 10 
ISTR{3) 21 
ISTR(4) 22 
ISTR(5) spillover 
Thru from previous 
ISTR(66) calculations 

Subsequent calls to calculational and report subroutines will 
transmit the information that only the first four elements of the ISTR 
array are to be used. The logic contained within subroutine ISTN has 
thus informed subroutine MATCLC (and subroutines HYDCAL, UNTREP and 
UNTENR which also utilize this information) in this example that the 
mass balance and related calculations are for the alkylation unit. This 
remains in force until the next call to subroutine ISTN is made. 
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The calls to subroutine ISTN are made from MATCLC and the main 
program. The calls from the main program are made to Identify the feed 
and product streams for the petroleum, shale and coal oil crude processing 
units. 

J. HYDCAL (AAA,SSS,HNN,PCTHY.NSTRM) 

This subroutine calculates the hydrogen content of process unit 

streams as a function of API gravity and sulfur and nitrogen content. 

The calling arguments are listed below: 

AAA specific gravity (converted in HYDCAL to API gravity) 

SSS weight percent sulfur content 

NNN weight percent nitrogen content 

PCTHY the weight percent hydrogen which Is returned 
to the calling subroutine 

NSTRM the number of active streams and thus the number 
of streams to be accessed within array ISTR. 

The calculations within HYDCAL are bypassed for any coded streams 
with an Identification number of 13 or lower since these are distinct 
light end hydrocarbon compounds (such as propane) whose hydrogen contents 
are known. The calculations are also bypassed If PCTHY has been prespecified 
and has a value other than the program Initialization value of 0.0, 

This allows model users the flexibility to make their own estimates of 
hydrogen content for the key aviation turbine fuel blending streams. 

The array STRMBP which is set within HYDCAL contains the mean 
boning points for each of the possible sixty-six process streams 
represented within the model (the first thirteen are set at 0.0 because 
these light end streams are not referenced). The BPI, APII and CKR 
arrays contain the correlation values obtained from the Technical Data 


Book - PetroleusEi Refining (1966), Figure 2B2.1, page 2-11, and refer to 
mean boiling point (®F), API gravity and carbon-to-hydrogen ratio. This 
correlation forms the basis of the HYDCAL prediction in which a double 
linear interpolation is performed between the stream API and mean boiling 
points, with reference to the tabular array values, to arrive at the 
estimated carbon-to-hydrogen ratio for each stream. Adjustments are 
then made for the sulfur and nitrogen contents to obtain the weight 
percent hydrogen content. A bias correction is programmed for shale and 
coal oil cuts to adjust the petroleum-based correlation which is used 
for these different hydrocarbon species. 

K. UNTREP {YIDARC,SULARC,NITARG,SPGARC,HYDARC,NSTR,NF, 

UNAME,UNAMEl ,UNAMEg)~ ^ 

This subroutine is called from MATCLC and reports the mass, sulfur, 
nitrogen and hydrogen balances for each refinery process unit. 

The report function is bypassed if the report level indicator IREP is 
set at less than 3. A set of preliminary calculations are always made 
in UNTREP,. however, since they involve the accumulation of mass, sulfur, 
nitrogen and hydrogen for the overall refinery material balance which is 
reported from the main program. These are accumulated in the arrays 
SVWGT, SVSUL, SVNIT, and SVHYD respectively, and are transmitted back to 
the main program thru the labeled common/OBAL/. 

The UNTREP subroutine arguments are listed below: 

YLDARC Stream masses 

SULARC Stream sulfur masses 

NITARC Stream nitrogen masses 

SP6ARC Stream speci'fic gravities 

HYDARC Stream hydrogen masses 

NSTR Total number of streams active in the ISTR array 

NF Total number of process unit feed streams 

UNAME The alpha numeric- 

UNAMEl process unit- 

UNAME2 name 
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The total mass and component masses are transmitted in units of 
pounds per day and converted to pounds per hour within UNTREP. Stream 
specific gravity values are converted to API gravity for reporting 
purposes. Subroutine UNTREP must distinguish between the process unit 
feed and product streams In order to report the material balances 
correctly. Since the CALL ISTN statement always lists feed streams 
first# the value NF, which provides the number of feed streams In the 
list# completes the information required by UNTREP. The use of the ISTR 
array containing the process unit stream code Identifiers permits the 
proper corresponding stream names to be referenced for report purposes. 

L. UNTENR (YLDARC,UNAM£#UMAME1#UNAME2) 

This Is a companion subroutine to the above which reports the 
refinery process unit energy balances if the report level Indicator IREP 
Is set at greater than E. A portion of the subroutine calculation 
proceeds 1n all cases and Involves the accumulation of energy related 
arrays for the overall energy consumption report Issued by the main 
program. These arrays are: 

SVBTU Btu*s consumed In the form of fuel# power, steam 
and hydrogen 

SVBTUS The above, converted to SI units 
SVFOE Expressed In barrels FOE 

SVPCT Energy consumption expressed as a percentage 

of total refinery Btu Input 

The UNTENR subroutine calling arguments are listed below: 

YLDARC The process unit masses. Only the first four# 
YLDARC(l) thru YLDARC (4) are referenced 

UNAME The alpha numerlc- 

UNAMEl process un1t- 

UNAME2 name 
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The first four stream identifiers in the YLDARC array always 

contain the following information for all process units: 

YLDARCnj Strem consumption 
YLDARC(2) Fuel consumption 
YLDARCf3) Power usage 
YLDARC(4) Hydrogen usage 

Steam, fuel and hydrogen usage are transmitted in units of barrels 
of fuel oil equivalent per day and converted to the reporting units 
within UNTENR. There is a potential for confusion between the hydrogen 
values for the catalytic reformer and hydrogen plant, which produce 
hydrogen, and all of the other process units which consume it. Therefore 
appropriate programming checks are made within subroutine UNTENR. 

M. BLENDF 

This subroutine provides the logic for calculating the sequence of 
fuel oil blends ranging from jet fuel aviation turbine blends to residual 
fuel oil. 

The first call from BLENDF is to subroutine BLDARY which sets up 

the volumetric and properties arrays for the entire forty-seven potential 

fuel oil blending components. The remaining BLENDF statements are 

involved with a series of call statements which produce each of the 

blends. A key to understanding the blending logic is to inspect the 

elanent of the VT(I) array to determine the blend being produced, where 

I is the blend code number. These code numbers are enumerated below; 

Blend 

Code 

Number Blend Description 

1 The total distillate fuel oil blend. 

2 The maximum distillate fuel oil blend that can be 
produced at SPECMD percent sulfur. 
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Bland 

Code 

Humber 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 
19 


Blend Description 

The remaining middle distillate fuel oil blend after 
blend 2 is produced from the total distillate fuel 
oil pool. 

The middle distillate blend of volu. VM]_ and percent 
sulfur SPMl. 

The middle distillate blend of volume VM2^ and percent 
sulfur SPH2 . 

The remaining middle distillate fuel oil blend after 
blends 4 and 5 have been produced. 

The total residual fuel oil blend. 

The maximum residual fuel oil blend that can be produced 
at SPECRF percent sulfur. 

The remaining residual fuel oil blend after blend 8 
is produced from the total residual fuel oil pool. 

The residual fuel oil blend of volume VRl and percent 
sulfur SPRl. 

The residual fuel oil blend of volume VR2 and percent 
sulfur SPR2. 

The remaining residual fuel oil blend after blends 10 
and 11 have been produced. 

The total middle distillate plus residual fuel oil 
blend after making all fixed volume blends, i.e. 
blends 4, 5, 10, and 11. 

The 525*F endpoint aviation turbine jet fuel blend. 

Blend 14 as limited by the minimum hydrogen content 
specification HDSPEC. 

The 650®F endpoint aviation turbine jet fuel blend. 

Blend 16 as limited by the minimum hydrogen content 
specification HDSPEC. 

The aviation turbine jet fuel blend at the specified 
variable endpoint EPSPEC . 

Blend IS as limited by the minimum hydrogen content 
specification HDSPEC. 
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The blend specification variables which are underlined above are 
defined on pagt'T 47 and 48 of Volume I. The general logic reflected in 
subroutine BLC " Is to produce all blends independently of one another 
wherever poss1b*i*e (this is not possible with blends 4,5, and 6 since 
6 Is defined as the blend remaining after blends 4 and 5 are produced, 
and similarly for blends 10, 11, and 12). 

A call is made to subroutine BLDREP from BLENDF in order to report 
the volume, composition and properties of each of the blends produced. 
These reports cannot be suppressed by user option. The program logic 
of producing and reporting all of the possible blend combinations yields 
the maximum information. The overall material balance reported from the 
main program includes only the foil owing: blend 19; the remaining 
middle distillate pool after making blend 19 (I.e. blend 1 minus blend 
19)r and blend 7. These are the aviation turbine jet fuel blends of 
specified endpoint and hydrogen content, the remaining middle distillate 
pool and the total residual fuel oil blend. These blends are of prime 
Interest in the overall material and volumetric balances and do not 
contain overlapping components, as do some of the other blends. 

A set sequence of subroutine calls Is made from within BLENDF for 
each of the fuel oil blends produced. 

These are: 

Subroutlne(s) 

SSORT or SSORTI 

BMAX 


Function 


Sorts blending components by maximuni sulfur 
or minimum hydrogen content. 

Calculates the maximum quantity of a blend 
of specified sulfur or hydrogen content that 
can be produced. 
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Sufarouti ne(s) 


Funct"! on 


BPROP Calculates the full set of blend properties. 

BLDREP Produces the blend reports. 

BMAX is omitted from the calling sequence for a particular blend if 
the particular blend volume is not limited by a sulfur or hydrogen 
specification (e.g. blends 1, 7, 14 and 16). 1 

N. BLDARY fIJK) 

BLDARY is called a total of three times from subroutine BLENDF. 

The first call establishes the stream volumes and all related properties 
for each of the potential forty-seven possible blending components. The 
second call adjusts these properties and stream volianes for the variable 
endpoint aviation turbine jet fuel blend (controlled thru the input 
£B?EC)„ and the third call restores- the original stream volumes for the 
middle distillate and residual fuel oil blending sequence contained 
within the final segments of subroutine BLENDF. The BLDARY argument, 

UK, is set at a value of 1 for the first and third call, and at 2 for 
the second call. 

The blend component volume and properties which are set within 
BLDARY are results of volumetric and mass calculations previously 
performed in the main program and subroutine MATCLC. The stream properties 
are established on the basis of internally stored data base property 
values contained in subroutines CBASE, CBASEl and DBASE. The following 
arrays are set up in subroutine BLDARY: 
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Array DescHption 


VOL 

Stream 

API 

Stream 

SUL 

Stream 

NIT 

Stream 

HDR 

Stream 

VIS 

Stream 

PAR 

Stream 

NAP 

Stream 

ARO 

Stream 

FRT 

Stream 

SPT 

Stream 

HTC 

Stream 

NUN 

Stream 

NST 

Stream 


volume, BPD 
specific gravity 
sulfur content, weight percent 
nitrogen content, weight percent 
hydrogen content, weight percent 
viscosity blending index 0 210®F 
paraffin content, weight percent 
napthene content, weight percent 
aromatics content, vfeight percent 
freezing point, ®F 
smoke point, millimeters 
heat of combustion, Btu per pound 
unit code (producing unit) 
name code 


The unit and name codes are established for the purpose of issuing 
blend reports and they are linked with a set of unit and stream alpha- 
numeric names which are referenced from subroutine BLDREP. Each element 


of the above arrays corresponds to a potential blending component, the 
identity of which is established in subroutine BLDARY. There are a 
total of forty-seven potential fuel oil blending components as listed 


below: 


Humber Fuel Blending Component Description 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
n 
12 

13 

14 

15 
15 

17 

18 

19 

20 


Of POOE 


purchased kerosene 

hydrotreated light kerosene - petroleum derived 
hydrotreated light kerosene - shale oil derived 
hydrotreated light kerosene - coal oil derived 
hydrocracker light kerosene - petroleum derived 
hydrocracker light kerosene - shale oil derived 
hydrocracker light kerosene - coal oil derived 
raw light kerosene - petroleum derived 
raw light kerosene - shale oil derived 
raw light kerosene - coal oil derived 
reformer endpoint control tower bottoms 
purchased heavy kerosene 

hydrotreated heavy kerosene - petroleum derived 
hydrotreated heavy kerosene - shale oil derived 
hydrotreated heavy kerosene - coal oil derived 
hydrocracker heavy kerosene - petroleum derived 
hydrocracker heavy kerosene - shale oil derived 
hydrocracker heavy kerosene - coal oil derived 
raw heavy kerosene - petroleum derived 
raw heavy kerosene - shale oil derived 
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Fuel Blending Component Description 


21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 


raw heavy kerosene - coal oil derived 

desulfurizer distillate from straight run gas oil 

desulfurizer distillate from/ coker gas oil 

desulfurizer light cycle oir‘^ 

desulfurizer light coker gas oil 

raw light cycle oil 

raw light coker oil 

purchased gas oil 

purchased residual fuel oil 

desulfurizer petroleum gas oil 

desulfurizer coker gas oil 

raw petroleum vacuum gas oil^ij 

thermal cracker fuel oil cut 

fluid cracker tower bottoms 

reduced crude - petroleum derived 

reduced crude - shale oil derived 

reduced crude - coal oil derived 

total coker gas oil 

heavy coker gas oil 

vacuum bottoms - petroleum derived 

vacuum bottoms ~ shale oil derived 

vacuum bottoms ~ coal oil derived 

visbreaker pitch 

raw shale gas oil 

raw coal gas oil 

hydrotreated shale gas oil 

hydrotreated coal gas oil 


The logic for adjusting the volume and properties of the jet fuel 
blending components comnences following statements 1314 and is contained 
within the DO LOOP for statement 10* The volume and properties are calculated 
by linear interpolation between those corresponding to 525 and 650®F, 


cut points which correspond to yields and properties given in the crude 
oil assays- The specified jet fuel blend endpoint, EPSPEC, is used to 
interpolate properties, 

0. S50RT (IL, lU) 

This subroutine sorts the volume and property arrays in order of 
increasing sulfur content for middle distillate and residual fuel oil 


(1) Petroleum derived. 
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blends. This sort is required in order to calculate the maximum quantities 
of fuel oil blends that can be produced with a specified sulfur content. 

The original component volumes and their corresponding properties are 
held in permanent arrays (for the duration of the run) which are set up 
in BLENDF, the subroutine which calls SSORT. These permanent arrays are 

necessary in order to reestablish the initial blend component values for 

\ 

subsequent blend calculations performed in BLENDF. The permanent blend 
arrays are identified below: 


Working 

Array 

Permanent 

Array 

Property 

VOL 

VHOLD 

Vol ume 

API 

AHOLD 

Specific gravity 

SUL 

SHOLD 

Sulfur 

NIT 

NIHOLD 

Nitrogen 

HDR 

HHOLD 

Hydrogen 

VIS 

VIHOLD 

Viscosity 

PAR 

PAHOLD 

Paraffins 

NAP 

NAHOLD 

Mapthenas 

ARO 

ARNOLD 

Aromatics 

SPT 

SPHOLD 

Smoke Point 

FRT 

FRHOLD 

Freezing Point 

HTC 

HCHOLD 

Heat of Combustion 

NST 

NSHOLO 

Stream Number 

NUN 

NUHOLO 

Producing Unit 


The subroutine arguments IL and lU define the range of blending 


components to be sorted and are selected from the blend component 
numbers defined in subroutine BLDARY. 

P. SSORTl (IL, lU) 


This subroutine is analogous to SSORT and differs only in that the 
sort is according to hydrogen content ranking as opposed to sulfur 
content. The working arrays are established in SSORTl in the order of 
decreasing hydrogen content. 
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Q, BMAX (lU, IL, IBL, IBT) 

This subroutine calculates the maximum volume of a blend of specified 
sulfur content, or hydrogen content, which can be produced from the fuel 
oil and jet fuel blend pools. 

The subroutine arguments IL and lU define the range of blending 
components to‘ be considered and are selected from the blend component 
numbers defined in subroutine BLDARY. Argument IBL is the blend number 
being calculated by BMAX (these blend identifiers are listed under 
the description for subroutine BLENDF). Argument IBT is the blend 
number from which the blend IBl is being produced. For example, the 
call statement preceding the DO 302 T66p in BLENDF is CALL BMAX (1, 11, 

19, 18) and conveys the following: 

0 Limit the calculation to VOL(l) thru VOL(ll) corresponding to 
the range from purchased kerosene to reformer feed preparation 
tower bottoms (this is the full range of jet fuel blend components 
as defined in subroutine BLDARY). 

0 BMAX will calculate the maximum quantity of aviation turbine 
fuel blend of hydrogen content HDSPEC - blend 19 as defined in 
subroutine BLENDF. 

0 The above blend is to be produced from blend 18 - the maximum 
jet fuel blend volume of endpoint EPSPEC. 

The composition of each blend calculated in subroutine BMAX is 
contained in array COMP (J, IBL), where J is the component number and 
IBL is the blend number, as defined in subroutines BLDARY and BLENDF, 
respectively. 
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The blend calculation will fail in BMAX if an impossible task is 
set forth, e.g. if all of the jet fuel blend components have a lower 
hydrogen content than HDSPEC. This situation triggers an error message 
and the run is then halted. Appropriate input adj’ustments must be made 
to obtain a successful run, such as lowering the blend specification or 
altering the refinery configuration in order to produce blending components 
of higher hydrogen contents. 

R. 8PR0P (JL, JU, ITl 

BPROP calculates the properties of each blend produced. The subroutine 
arguments denote the range of blend components considered for a particular 
blend (JL and JU) as defined in subroutine BLDARY, and the blend number 
(IT) as defined in subroutine BLENDF. The following arrays appearing in 
BPROP contain the blended properties of each blend. The method of 
blending is indicated: 


Array 

Description 

Blending Method 

AP (IT) 

Specific gravity 

volumetrically 

SU (IT 

Sulfur content 

by weight 

NI (IT) 

Nitrogen content 

by weight 

HD - (IT) 

Hydrogen content 

by weight 

PA (IT 

paraffin content 

by weight 

flA (IT 

napthene content 

by weight 

AR (IT) 

aromatic content 

by weight 

HT (IT 

heat of combustion 

by weight 

VI (IT 

viscosity blending index 

vol umetri cal ly 

SP (IT) 

smoke point 

by reciprocal 

volumetric 

blending 

FR (IT) 

freezing point 

by blending 
function BFPF 
as defined in 
BPROP 


Subroutine BPROP is called from BLENDF. 
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S, BLDREP (IL, lU, JB) 

This subroutine is called from BLENDF in order to report each of 


the blends produced. The blending component names are referenced thru 
the array SNAHES as defined in the BLOCK DATA subroutine. The refinery 
process units which produce these jet fuel* distillate and residual fuel 
oil blends are also reported using the UNAMES array. Speci,f1c gravity 

I 

is converted to API gravity and both are reported. The blending reports 
issued by subroutine BLDREP are reported for all values of the report 
option indicator IREP. The BLDREP subroutine arguments denote the range 
of blend components considered for a particular blend (IL and III) and 
the blend number (JB). 

T. CONAPI (JJJ) 


This subroutine is called to convert all streams (API) gravities to 
specific gravities (JJJ=1) for the calculation of the process unit 
material balances and from specific gravity to API gravity (JJJ=2) in 
order to restore the initial values for reporting purposes. 

U. ECON 


This subroutine calculates the 
report option indicator IREP is set 
order of calculation is as follows: 
Description 

Process unit Investments 
Nelson refinery complexity factor 
Refinery offsites investment 
Investment report 
Refinery fixed '^operating cost 
calculation and reporting 
Refinery variable operating 

cost calculation and reporting 


overall refinery economics if the 
at greater than a value of 1. The 

Subroutine Location 

preceding statement 129 
statement 129 thru 14 
statement 14 thru 148 
statement 148 thru 48 

\ 

statement 48 thru 4 
statement 4 thru 8 
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Description 

Crude oil cost calculation and 
reporting 

Refinery product revenue calculation 
and reporting 

Summary refinery economic calculations 
and reporting 


Subroutine Location 


statement 8 thru 12 
statement 12 thru 529 
statement 529 thru 801 


The above statement numbers refer to the ECON subroutine labels. 


which are not numbered sequentially. 


2.2. Subroutine Heirarchy 

The heirarchy of program subroutines is described below. The 

underlined subroutine calls those which are listed below it, 

WAIN PROGRAM 

BLOCK DATA 

INITB 

INITV 

UBASE 

CBASE 

CONAPI 

ISTN 

MATCLC 

BLENDF 

ECON 

MATCLC 

ISTN 

HYDCAL 

tlNTREP 

UNTENR 

BLENDF 

BLDARY 

SSORT 

SSORTI 

BPROP 

BMAX 

BLDREP 


CBASE 

CBASEl 
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3,0 COMMON BLOCK DESCRIPTIONS 


3.1 Common Blocks Crossreferenced by Subroutine 


Common Block/ 


BLEND CB HTCOM ISTNUH UNDEN IB 


Subroutfne 


REFMOD (Main) 

ECON 

CBASE 

CBASEl 

HYDCAL 

BLOCK DATA 

CONAPI 

ISTN 

UBASE 

INIT8 

UNTREP 

UNTENR 

BLDARY 

BLENDF 

SSORTI 

BLDREP 

BPROP 

SSORT 

BMAX 

HATCLC 

INITV 


XXX 

X 

X 

X 

X X 

X X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X X 



X X 

X 
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xxxxxxxx 


Common Block/ INAME MATCAL TITLER ML lii M JH. M UB6 


Subroutine 


REFHOD (Main) 

ECON 

CEASE 

CBASEl 

HYDCAl 

BLOCK DATA 

CONAPI 

ISTN 

UBASE 

INITB 

UNTREP 

UNTENR 

BLDARY 

BLENDE 

SSORTV 

BLDREP 

BPROD 

SSORT 

BKAX 

HATCLC 

INITV 


X X 

X 


X X X X X X X 

X 



X 


X 


X X X X X 

X X X X 


X X 


X X 


X !x 
X 


X X X x X X 


X 


X 

X 


X 

X 


X 

X 


X 

X 


X 


X 


ComKon Block/ 


PASER EBAL HYDROS OBAL INITB 


ECONOM 


Subroutine 

REFMOD (Main) 

ECON 

CBASE 

CBASEl 

HYDCAL 

BLOCK DATA 

CONAPI 

ISTN 

UBASE 

INITB 

UNTREP 

UNTENR 

BLDARY 

BLENDE 

SSORT.1 ■ 

BLDREP 

SSORT 

BHAX 

MATCLC 

INITV 


X 

X 


X 

X 

X 


X 

X 


X X 
X 


X 

X 


X 

X 


X 


3.2 Description of Labeled Common Blocks 


BLEND 


HTCOM 


ISTNUM 


LENDEN 


IB 


INANE 


HATCAL 


TITLER - 

UBl, UB2, - 
UB3 


This coiMon block contains all of the blend property 
arrays for the jet fuel, distillate and residual fuel 
oil blends. This block also contains the blend volumes, 
compositions and the number of blend components. 

This common block contains all of crude assay properties, 
the specified crude volumes, and the total number of 
crudes in the data base. 

Btu conversion factors are contained in this common 
block. HCON(l) = Btu per lb of steam; HC0N(2) = Btu 
per lb fuel gas; HC0N(3) =* Btu/kWh; HC0N(4) = Btu per 
lb hydrogen. This block also contains the calculation of the 
Btu per day content of the total refinery oil feed. 

This coninon block contains the array of active internal 
stream numbers- pertinent to a given calculation and also 
the total number of these streams. 

This common block contains the specific gravity of the 
first thirteen streams which correspond to distinct 
hydrocarbon species such as propane and propylene. 

This common block contains unit capacities, stream 
dispositions between alternate units, process unit 
severities, blending stock purchases, byproduct sales 
volumes and blend specifications ~ all of which are 
initialized in subroutine INITB. 

This common block contains the array of crude oil 
names, product stream names and process unit names 
- all of which are used for reporting purposes. 

This common block contains the volumetric (BPD) 
and weight flow (pound per hour) rates of all feed 
and intermediate process unit streams. 

This common block contains the current case run alph- 
numeric title. 

These common blocks contain the physical properties 
of process unit feed and product streams. These 
include specific gravity; sulfur, nitrogen and 
hydrogen contents; and viscosities for a portion of 
the process units. Also included are the process 
unit fractional yields. 
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UB4 

UBS 

UB6 

PAGER 

EBAL 

HYDROS 

OBAL 

INITB 

ECONOM 


This common block contains the remaining stream vis- 
cosities; intermediate stream clear and leaded octane 
values and the differential specific gravities between 
process unit feed and product streams. 

This common block contains the limiting catalytic cracker 
gas oil conversion values and the limiting catalytic 
cracker platformate clear octane values which mark the 
severity range. 

This common block contains the special-light and heavy 
kerosene properties which are pertinent to jet fuel 
blending. These include heat of combustion, PNA*, 
smoke and freeze points. This common block also contains 
the differential PNA values between process units feed and 
product streams as a function of hydrotreating severity. 

This common block contains the page number Indicator 
and the report indicator IREP which controls the level 
of detail which is reported. 

This common block contains the oil refinery energy 
usage accumulators for steam, fuel, power and hydrogen. 

This common block contains the calculated hydrogen 
contents for all the process unit internal streams. 

The hydrogen content is calculated as a function of 
specific gravity and mean boiling point, unless otherwise 
specified. 

This common block contains accumulators for the sulfur, 
nitrogen and hydrogen contents of all oil refinery 
feed and Intermediate streams, as well as the mass 
accumulator. These are used to calculate the refinery 
mass, sulfur, nitrogen and hydrogen balances. 

This common block contains oil refinery process unit 
feed and key product volumes. These are initialized 
In subroutine INITB before each case run calculation. 

This common block contains all process unit fixed 
and variable costs; feed and product values; the refinery 
stream factor; scale factors used in the calculation 
of process unit investment costs; the capital investment 
carrying charge; and the cost of electricity. 


PNA refers to the weight percent of paraffins, naphthenes, and aromatics 
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4.0 FUNCTION STATEMENT DESCRIPTIONS 


MAIN PROGRAM 


BFPF (FREPT) 

The conversion from freezing point (®F) 
to blending index value 

FREP (BIFRP) 

The conversion from blending index value 
to freezing point (*F) 

CSl(X) 

Conversion from 42 gallon barrels to 
cubic meters 

CS2 (X) 

Conversion from pound per hour to 
kilogram per second 

CS3 (X) 

Conversion from cubic feet to cubic 
meters 


B. SUBROUTINE BIDARY 

BFPF (FREPT) The conversion from freezing 

to blending index value 

FREP (BIFRP) The conversion from blending 

valu»^: to frspzing point (®F) 

TAPI (X) The cvrtv»:r.’‘S!on from specific 

API 

C. SUBROUTINE BPROP 

BFPF (FREPT) See above 

FREP (BIFRP) See above 

D. SUBROUTINE HATCLC 

TAPI (XXX) Conversion from specific gravity 

in API gravity ' 

TSPSG (yyy) Conversion from API gravity to 

specific gravity 


point (‘*F) 
index 
gravity to 
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5.0 PROGRAM VARIABLE DESCRIPTIONS 


The following sections provide descriptions of the symbolic variables 
appearing in labeled common blocks and local data statements. The 
elements contained in the dimensioned stream related arrays, such as 
STR, are in the order given in Exhibit 5-Internal Refinery Stream and 
Process Identifiers -page 73 of Volume I (NASA CR-135333). Jet fuel, 
distillate fuel oil and residual fuel oil blend component arrays, such 
as VOL and API, contain elements in the order of the fuel blending 
components description for the 47 possible components as given in 
Section 2.1 of this volume under the description of subroutine BLDARY. 

Fuel blend arrays relating to the total jet fuel and fuel oil blends, 
such as VT and AP, contain elements In the order given under the Section 

2.1 description of BLENDF, which enumerates 19 possible fuel blends in 
the sequence in which they are entered into the arrays. 

5.1 Input Variables 

AIT input variables are fully described in Volume I (NASA CR- 
135333) on pages 22 thru 54. Further description would be redundant. 

The input data variable names described are the same as the program 
symbolic names. 

5.2 Additional Variables Appearing in Common 

The following common variables are not input variables and are 
described below under their respective labeled common block. The variable 
type is indicated: R is real; I Is integer; RA is real array; and lA is 
integer array variable. 
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A. CQMHON/BLEND/ 


VARIABLE 

TYPE 

DESCRIPTION 

VOL 

RA 

Array containing volume (BPO) of each 
jet fuels middle distillate and residual 
fuel oil blending component. 

API 

RA 

Corresponding specific gravity. 

NIT 

RA 

Corresponding weight percent nitrogen 
content. 

SUL 

RA 

Corresponding weight percent sulfur 
content. 

VIS 

RA 

Corresponding viscosity blending index 
value @ 210®F. 

HDR 

RA 

Corresponding hydrogen content in weight 
percent. 

FRT 

RA 

Corresponding freezing point in degrees F. 

PAR 

RA 

Corresponding paraffin content, weight 
percent. 

NAP 

RA 

Corresponding naphthene content in weight 
percent. 

ARO 

RA 

Corresponding aromatic content in weight 
percent. 

SPT 

RA 

Corresponding smoke point in millimeters. 

HTC 

RA 

Corresponding heat of combustion in 
BTU per pound. 

NST 

lA 

Corresponding stream name identifying 
reference number. 

NUN 

lA 

Corresponding identifying reference 
number for the process unit producing 
the stream. 

VT 

RA 

The array containing the total volume 
of each of the possible 19 blends. 
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VARIABLE 

TYPE 

DESCRIPTION 

AP 

RA 

The corresponding specific gravity of 
each of the blends. 

NI 

RA 

Corresponding weight percent nitrogen 
content. 

SU 

RA 

Corresponding weight percent sulfur 
content. 

VI 

RA 

Corresponding viscosity blending index 
at 210®F. 

HD 

RA 

Corresponding weight percent hydrogen 
content. 

FR 

RA 

Corresponding blend freezing point in 
degree F. 

PA 

RA 

Corresponding weight percent paraffins 
of each blend. 

NA 

RA 

Corresponding weight percent naphthenes 
of each blend. 

AR 

RA 

Corresponding weight percent aromatics 
of each blend. 

SP 

RA 

Corresponding smoke point of each blend 
in millimeters. 

HT 

RA 

Corresponding heat of combustion of 
each blend in BTU per pound. 

PCT 

RA 

Total percentage volume of each blend. 
This will equal 100.0 with a negligible 
round-off discrepancy. 

COMP (1,0) 

RA 

The volume percent of the I‘th component 
in the O'th blend. 

NBC 

I 

The total number of blend components. 
Set at 47, 

B. COMMON/CBA 



VARIABLE 

TYPE 

DESCRIPTION 

NC 

I 

Total number of potential data base 
crudes, set at 35. 
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C> COMMQN/HTCQM/ 


VARIABLE 

npE 

DESCRIPTION 

HCON 

RA 

The heat content of steam* fuel gas, 
electricity and hydrogen, respectively. 
Expressed in BTU per pound, except for 
electricity in BTU/kWH, 

BTUPOT 

t-- 

R 

Total refinery heat content input in 
BTU per day. 

D. COMMON/ ISTNUM/ 

t 

VARIABLE 

TYPE 

DESCRIPTION 

T?TR 

lA 

The array of stream identification code 
numbers which are currently active in 
the program calculation. 

NSTOT 

I 

The total number of active streams within 
array ISTR. 

E. COMMON/LENDEN/ 


VARIABLE 

TYPE 

DESCRIPTION 

SPGLE 

RA 

The array containing the specific gravity 
of light and components. 

F. COMMON/ INAME/ 


VARIABLE 

TYPE 

DESCRIPTION 

CN 

RA 

The array containing the data base crude 
oil names. 

SKAMES 

RA 

The array containing the process unit 
stream names. 

UNAMES 

RA 

The array containing the process unit 
names. 

G. COMMON/MATCAL/ 


VARIABLE 

TYPE 

DESCRIPTION 

SIR 

RA 

The array containing the BPD volume of 
each refining stream. 

STC 

RA 

The array containing the volume of petroleum 
based crude oil streams. 

SIRS 

RA 

The arrav conl!ainino the volunt^ of oil 

crude unit streams. 
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STRC 

RA 

The array containing the volume of coal 
o1l crude unit streams. 

STHK 

RA 

The array for kerosene hydrotreater 
feed and product streams. 

STCO 

RA 

The array for fluid coker feed and 
product streams. 

STVB 

RA 

The visbreaker stream array. 

STCCG 

RA 

The fluid catalytic cracker stream 
array for virgin gas oil feed. 

STCCC 

RA 

The fluid catalytic cracker stream 
array for coker gas oil feed. 

STCCX 

RA 

The fluid catalytic cracker stream 
array for total combined feed. 

sncc 

RA 

The thernal cracker stream array. 

STHYGG 

RA 

The gas oil hydrocracker stream array 
for virgin gas oil feed. 

STHYGC 

RA 

The gas oil hydrocracker stream array 
for coker gas oil feed. 

STHYGX 

RA 

The gas oil hydrocracker stream array 
for toal combined feed. 

STGDG 

RA 

The gas oil desulfurizer stream array 
for virgin gas oil feed. 

STGDC 

RA 

The gas oil desulfurizer stream array 
for coker gas oil feed. 

STGDX 

RA 

The gas oil desulfurizer stream array 
for total combined feed. 

STDDL 

RA 

The distillate desulfurizer stream array 
for light cycle oil feed. 

STOOD 

RA 

The distillate desulfurizer stream array 
for virgin distillate desulfurizer. 

STDDX 

RA 

The distillate desulfurizer stream array 
for total combined feed. 
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STDHD 

RA 

The distillate hydrocracker stream array 
for distillate feed. 

STDHL 

RA 

The distillate hydrocracker stream array 
for light cycle oil feed. 

STDHX 

RA 

The distillate hydrocracker stream array 
for total condsined feed. 

STREP 

RA 

The catalytic reformer stream array. 

STBI 

RA 

The butane isomerizer stream array. 

STAL 

RA 

The alkylation unit stream array. 

STPO 

RA 

The polymerization unit stream array. 

STH2 

RA 

The hydrogen plant stream array. 

STSia.1 

RA 

The shale oil hydrotreater stream array 
for light kerosene range feed. 

STSKHl 

RA 

The shale oil hydrotreater stream array for 
heavy kerosene feed. 

STSKXl 

RA 

The shale oil hydrotreater stream array 
for combined feed. 

STCKLl 

RA 

The coal oil hydrotreater stream array for 
light kerosene range feed. 

STCKHl 

RA 

The coal oil hydrotreater stream array 
for heavy kerosene feed. 

STCKXT 

RA 

The coal oil hydrotreater stream array 
for combined feed. 

STSGHl 

RA 

The shale oil gas oil hydrotreater stream 
array. 

STCGHl 

RA 

The coal oil gas oil hydrotreater stream 
array. 


All of the above MATCAL common variables have companion variables 
beginning with W rather than S, These are the corresponding streams 
expressed in mass flow (pounds per hour) as opposed to volume flow 
(barrels par day). 
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'i 


H, COMHON/ECONOM/ 


VARIABLE 

nPE 

DESCRIPTION 

DIFF 

R 

The net refinery fuel which must be 
purchased In BFOE per day. 

I. COMHON/PAGER/ 


VARIABLE 

TYPE 

DESCRIPTION 

IPAGE 

I 

The current report page number. 

J. COMMON/EBAL/ 


VARIABLE 

TYPE 

DESCRIPTION 

SVBTU 

RA 

The total refinery energy consumption 
of steam, fuel, power and hydrogen, 
respectively, expressed In BTU per 
hour. 

SVBTUS 

RA 

The above in joules per second. 

SVFOE 

RA 

The above In barrels FOE per hour. 

SVPCT 

RA 

The above as a percentage of total 
refinery heat content Input. 

K. COMMON/HYDROG/ 


VARIABLE 

TYPE 

DESCRIPTION 

HYDVAL 

RA 

The hydrogen content In weight percent 
of each refinery stream. 

L. COHMON/OBAL/ 



VARIABLE TYPE DESCRIPTION 


SVSUL 

RA 

The total refinery sulfur by stream In 
pounds per hour. 

SVNIT 

RA 

The total refinery nitrogen by stream In 
pounds per hour. 

SVHYD 

RA 

The total refinery hydrogen by stream 
In pounds per hour. 

SVWGT 

RA 

The total mass by stream 1n pounds per hour 
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5.3 Local Variables Appearing in Data Statements 

The key local variables appearing In data statanents are described 
below and the subroutine In which they are defined is Identified. 


SUBROUTINE 

VARIABLE 

TYPE 

DESCRIPTION 

Main Program 

CONI 

R 

Conversion factor, pounds of water/ 
42 gal barrel O 60F«349.788. 

ECON 

RPLEX 

RA 

Nelson Complexity points on graph of 
complexity vs. offsites cost. 

ECON 

OFFSP 

RA 

Nelson offsite costs expressed as a 
percentage of total refinery Investment. 

ECON 

CONI 

R 

Conversion factor, cubic meters/42 
gallon barrel = 0.158987. 

ECON 

mz 

R 

Cubic meters/cublc foot = 0.02832, 

ECON 

C0N3 

R 

Gallons/cublc foot = 7. 48052- 

ECON 

C0N4 

R 

Barrel s/cubic foot = 0.1781076 

ECON 

GALS 

R 

Conversion factor, gallons/barrel=42.0 

ECON 

UNTXWF 

RA 

Nelson Con^lexlty factors for each 
refinery process unit In order of 
Exhibit 5, p. 75, Volume I. 

HYDCAL 

BPI 

RA 

Array of API Technical Data Book mean 
boiling points taken from nomograph 
correlation of carbon- to-hydrogen 
ratio as a function of API and mean 
boiling points. 

HYDCAL 

BLDMBP 

RA 

Array of jet fuel and fuel o1l mean 
boiling points (In order of VT array). 

HYDCAL 

STRMBP 

RA 

Array of fuel oil blending component 
mean boiling points (1n order of VOL 
array) . 

HYDCAL 

HYDLE 

RA 

Array of light end components carbon- 
to-hydrogen ratios. 
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SUBROUTINE 

VARIABLE 

TYPE 

DESCRIPTION 

HYDCAL 

APII 

RA 

Array of API gravities corresponding 
to the BP I array. 

HYDCAL 

CHR 

R 

Array of carbon-to-hydrogen ratios 
corresponding to the BPI array. 

UNTENR 

BFOE 

R 

Conversion constant, BTU per barrel 
FOE = 6,05 million. 

UNTENR 

CONI 

R 

Pounds of water per 42 gallon barrel 
@ 60"F = 349.788. 

UNTENR 

C0N2 

R 

Kg/second per pound/hour = 1.25997 E-04. 

UNTENR 

C0N3 

R 

Ooules/kg per BTU/pound - 2326.0. 

UNTENR 

C0N4 

R 

Joules/second per BTU/hour = 0.29307, 

UNTENR 

CONS 

R 

Joules per BTU @ 60‘’F ^ 1054.68, 

BLDARY 

BLDHBP 

RA 

Same as BLDMBP in HYDCAL. 

BLDREP 

CONI 

R 

Cubic meters per 42 gallon barrel = 
0,1589873. 

BLDREP 

C0N2 

R 

Density of water @ 60®F = 999.04 
kg/cubic meter. 

BLDREP 

CONS 

R 

-Joules/ kg per BTU/pound =2326.0. 

MATCLC 

CONI 

R 

Pounds water per barrel @ 60®F = 
349.788. 

9 

UNTREP 

CONI 

R 

Pounds of water per 42 gallon 
barrel = 349.788, 

UNTREP 

CQN2 

R 

Kg/second per pound/hour * 
1.25997E-04. 

UNTREP 

CONS 

R 

Density of water, at 60®F * 999.04 
Kg per cubic meter. 


6.0 DESCRIPTIOH OF PROGRAM FILES 


The data bases for crude oil assay data, process unit yields, 
stream properties and economic data are self-contained within the Computer 
Model For Refinery Operations. Therefore, the program is not interfaced 
with external data base files and the use of files is limited to the 


standard input/output files and small temporary disk files which are 

These files are described below: 


purged after run termination- 

Logical 
• File 

Number Type 

5 Card Reader 

£ Printer 

17 Disk 

18 Disk 

16 Disk 

19 Disk 


( 1 ) 

Function 

Program input stream 

Printed program output 

Store user directed crude oil 
assay data base changes for 
the duration- of the run (entered 
under reader card &CDATA)- 

Store process unit data base 
changes entered under header 
card &UDATA. 

Store process unit data base 
changes entered under header 
card &UDATA1. 

Store card images for the current 
run problem entered under header 
card &PDATA. 


20 


Disk 


Store the above information 
on File 20 if the run problem 
corresponds to a base case. 


(1) All files are accessed from the main program only, with the exception 
of File 6 which is referenced from the main program and from the 
report subroutines. 
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Files 16, 17, and 18 are reread In the main program In order to 
overlay data base values stored within subroutines CBASE, CBASEl and 
UBASE. Files 19 and 20 are used in tandem. The base case problem Input 
Is stored on File 20 and the subject problem changes are stored on File 
19. File 19 Is reread following a read of file 20. Therefore, base 
case changes are not additive from subject case to subject case, but are 
always made with reference to the starting base case of a case study 
sequence. For a stand alone case, only file 19 Is used. 
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7.0 PROGRAM SIZE LIMITATIONS 


Tha Computer Model For Refinery Operations requires approximately 
400K bytes of storage on the IBM 360/67 computer and about 225 K words of 
storage on the CDC 6600 computer. Essentially all of the dimensioned 
array space reserved In the program is utilized. The key array limitations 


are outlined below: 

ITEM MAXIMUM 

Number of refinery process streams 66 

Number of refinery process units 20 

Number of crude oils (all types) 35 

Number of jet fuel and fuel oil blending 

components (total for all blends) 50 

Number of distinct jet fuel plus distillate 

and residual fuel oil blends 20 

Total number of items of crude oil assay 
information (individual cut fraction 
yields, properties, etc,) 57 

Number of hydrotreater severity levels for 

each shale and coal oil process units 3 

Number of purchased fuel oil blending stocks 4 
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8.0 CONCLUDING REMARKS 


This report has presented detailed systems and programming documentation 
for The Computer Model For Refinery Operations. The report is Volume 
III of three volumes covering the description, application, and documentation 
of the refinery calculation program. This volume is intended to be used 
in conjuncipion with Volume I and II and with the program Fortran listing. 

The information represents a detailed description of all program subroutines, 
common blocks, program function statements, data statements and program 
files used during run execution. The preceding volumes (NASA CR-1 35333) 
and NASA CR-1 35334) are available from the NASA Project Manager. Computer 
tapes can be purchased through the Computer Software and Management 
Information Office (COSMIC), 112 Barrow Hall, University of Georgia, 

Athens, QA 30602, under the number LEW-13047. 
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